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Abstract: We consider Next-to-Leading-Order QCD corrections to ADD graviton ex- 
change relevant for Deep Inelastic Scattering experiments. We calculate the relevant NLO 
structure functions by calculating the virtual and real corrections for a set of graviton in- 
teraction diagrams, demonstrating the expected cancellation of the UV and IR divergences. 
We compare the NLO and LO results at the centre-of-mass energy relevant to HERA ex- 
periments as well as for the proposed higher energy lepton-proton collider, LHeC, which 
has a higher fundamental scale reach. 



Contents 



i 


Introduction 


i± 




Leading order calculation 




12. 11 Structure functions definition 


i 




2.2 Sum rules 





i 


Graviton exchange NLO results 






STTI Quark constituent contribution 


I 




3^21 Gluon constituent contribution 


i 


i 


Interference with SM gauge bosons 


12 


i 


Numerical results 


17 




Conclusions 


20 



1. Introduction 

Models with extra dimensions have been proposed to address the hierarchy problem. These 
can in general be classified into two groups depending on the geometry of the extra di- 
mensions; flat or factorisable models and warped or non-factorisable models. The former 
includes the model proposed in ^ by Arkani-Hamed, Dimopoulos and Dvali (ADD) and 
its variants, while the latter includes the Randall-Sundrum (RS) model |p. 

The large extra dimensions model introduced in |||, ^ suggests that the Standard 
Model (SM) particles live in the usual S+l-dimensional space, while gravity propagates 
in a higher D-dimensional space. The weakness of gravity with respect to other forces is 
imposed by the large size of the compactified extra dimensions. In this model Newton's 
constant is expressed as 

= mI = 87rii'^M|+", (1.1) 

where Mp is the Planck mass. Mo ~TeV is the fundamental mass scale, n is the number 
of extra dimensions and R the radius of the compactified space, assumed to be a torus in 
this model. The size of the extra dimensions is determined by their number, n, and the 



fundamental mass scale, M^, through Eq. 1.1. For a fundamental scale of the order of 



1 TeV and one extra dimension, its size would be of the order of the solar system size, 
causing large deviations from the inverse square law of gravitation. Therefore one extra 
dimension is already ruled out by experiments. For more extra dimensions, the size gets 
rapidly smaller, not contradicting submillimetric gravity measurements Q. 
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In the ADD model Q , the graviton corresponds to the excitations of the D-dimensional 
metric. These can be expressed as a tower of Kaluza-Klein (KK) modes. The interaction 
Lagrangian for gravitons and the SM fields is given by 



1 



£,„t = -^VGlT^^ (1.2) 



with Mp = Mp/ vSvr the reduced four-dimensional Planck mass, T^jy the energy momen- 
tum tensor of the SM fields and n = (ni, n2, n„) a n-dimensional vector of integers 
labelling the massive gravitons. A set of Feynman rules for the interactions of these mas- 
sive KK states with the Standard Model particles was presented in |6|. Note that this 
is an effective theory of interactions valid in the cisplanckian region, that is when the 
centre-of-mass energy of the parton collision is smaller than the fundamental scale Mp). 
The predictions of this effective theory fail as we approach the quantum gravity scale. As 
we reach the planckian region where y/s ~ Md, we need theoretical input from quantum 
gravity, and so it is presently impossible to predict experimental signals reliably. 

The KK resonances have masses equal to m(^) = \n\/R. This results in a small mass 
gap of order R~^, e.g. for one extra dimension of size 1 /xm the mass gap is 0{1 eV). This 
renders different masses practically indistinguishable and permits replacing the sum over 
discrete mass values by an integral over a continuum with a given density of states. This 
density of states is obtained by considering the number of modes with KK index between 
\n\ and |n| -|- dn, and is given by 

dN Mp s~i n ox 

P^'^^ = d^ = ^'''J^"' ■ ^'-^^ 

As seen from the Lagrangian in Eq. graviton interactions are suppressed by inverse 
powers of Mp. Nevertheless, the summation over the large number of accessible KK modes, 



that is the integration over the density in Eq. |l.3| , cancels the dependence on Mp. As an 



example, the inclusive graviton production cross section, am, is expected from the graviton 
couplings to be proportional to Mp , which combined with Eq. |1.3| exactly cancels the 
dependence on Mp. This leads to an effective interaction suppressed by inverse powers of 
the fundamental mass scale Mp), thus giving observable effects for Mpj near the TeV scale. 

Interesting phenomenological implications of this extra dimensions model have been 
studied extensively in the literature. These include real graviton emission and virtual 
graviton exchange. A set of processes, such as graviton plus gauge boson production, was 
studied in [^, |^, where the Feynman rules were first presented. The experimental signature 
for graviton production is missing energy, as decay into SM particles is suppressed by a 
factor of 1/Mp, which is not compensated by phase space. Therefore gravitons behave 
as heavy and stable particles, once produced. One should expect missing energy signals 
at both lepton and hadron colliders. This missing-energy signal does not correspond to a 
fixed invisible-particle mass as the graviton has a continuous distribution in mass. This 
differentiates graviton searches from other Beyond the Standard Model (BSM) physics 
searches, such as supersymmetry which would also give a missing energy signal. The 
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experimental signatures depend strongly on the number of extra dimensions, n, and the 
fundamental mass scale, Mp, and so we should be able to determine or at least constrain 
both of these model parameters. Virtual graviton effects will cause deviations from SM 
predictions for fermion and boson pair production. The cross sections are divergent at tree 
level, and thus one has to introduce an ultra-violet cut-off that is usually taken to be of the 
order of the fundamental scale Q . The search for virtual graviton effects is complementary 
to the production search, and will shed light on this cut-off and its dependence on the 
fundamental scale and the number of extra dimensions. 

Both the LEP and Tevatron experiments have set limits on Md. The combined results 
of the LEP experimental collaborations for the monophotons channel are summarised in Q. 
The CDF and DO collaborations have searched for signals in the monojets and monophotons 
channels for real emission Q as well as dilepton and diphoton channels for virtual exchange 
effects. For a brief summary of the experimental searches and the current constraints on 
the fundamental scale and the number of extra dimensions see Q and references therein. 
A more recent study of the dijet angular distribution in DO was performed in |1C]. The 
most stringent current constraints on the fundamental scale are set at 1-1.5 TeV, but these 
also depend on the number of extra dimensions. The LHC is expected to probe extra 
dimensions effects up to higher scales, due to the larger centre-of-mass energy. Recently, a 
first study of the bounds set on the cut-off scale from the dijet cross section measured at 
the LHC H, m was presented in which are in the few TeV region. 

Recently, several graviton processes were calculated at Next-to-Leading Order (NLO) 
in QCD. Leading order calculations suffer from uncertainties due to the choice of renor- 
malisation and factorisation scales. NLO results are less dependent on the scale choice 
and thus improve theoretical predictions. Calculations have been performed for both ADD 
and RS models. Here we just mention a few of the processes studied at NLO QCD: gravi- 
ton production [0, |l^, graviton and photon associated production [jl6|, graviton plus two 
jets production |17], heavy resonance graviton production and decay into top quarks |18|, 
diphoton production ||l^, graviton plus Z production pO|, Drell-Yan pit] and gravitational 



scattering at transplanckian energies [22|. 

Graviton contributions to Deep Inelastic Scattering (DIS) have been considered by 
HERA experiments. The absence of any deviation from the SM prediction was used to set 
limits on the fundamental mass scale. The most recent limits for e~^p collisions from the 
HI experiment are given in [23|, while more information about the analysis method can be 
found in |24]. These studies use the leading-order (LO) prediction for the signal and NLO 
QCD corrections are not taken into account. Given the current limits on the fundamen- 
tal scale and the centre-of-mass energy available at the HERA experiments, discovery of 
graviton effects does not fall within the experimental reach. 

The graviton contribution to DIS will be more important in the context of the proposed 
Large Hadron electron Collider (LHeC) [^5|. The proposed design suggests using the LHC 
ring and a linear electron accelerator to collide 7 TeV protons with 140 GeV electrons. 
The significant increase of the centre-of-mass energy will allow much higher momentum 
transfers to be probed, with the effect of graviton DIS becoming important and much 
higher mass scales becoming accessible. 
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Figure 1: Diagrams for q£ and gt — gt. Red denotes the lepton line. 



In this paper, we evaluate the real and virtual corrections for simple graviton exchange 
relevant for DIS experiments. We present analytic results for the relevant structure func- 
tions at LO and NLO. We then present numerical results on the differential cross section 
and comment on the effect of NLO QCD corrections at HERA and LHeC energies. 



2. Leading order calculation 
2.1 Structure functions definition 

The leading order results for graviton exchange in DIS experiments were first published 



in [26 1 soon after the Feynman rules for gravitons became available. The corresponding 
parton-level Feynman diagrams for graviton exchange are shown in Fig. |l|. We note that 
there is a LO gluon diagram contributing to the scattering, in contrast to photon mediated 
DIS where gluon constituent scattering can only occur at NLO QCD. 

We start by considering the partonic differential cross section for the exchange of a 
graviton between a quark and a lepton, given by: 

(fa vr Q4 



dQ2 32M| y2 



[32 - 64y + 42y^-102/^ + /], (2.1) 



with and y the usual DIS variables and Ms the cut-off introduced to regularise the UV 
divergence of the graviton exchange. To obtain the hadronic cross section, the partonic 
cross section is convoluted with the proton parton distribution functions which in this case 
are sampled at momentum fraction x, giving the double differential cross section: 



4 



where the sum is over all quark flavours, with each quark flavour contributing equally. For 
comparison purposes we write down the corresponding expression for photon exchange in 
unpolarised deep inelastic scattering using the same variables: 

.[y^ -2y + 2]y-el[q{x) + q{x)i (2.3) 



where the contribution of each flavour depends now on the electric charge squared. The 
difference between the two expressions originates from the different dimension operators 
describing the interactions of photons and gravitons with the quarks. It is customary to 
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re-express the cross section for photon exchange for electron DIS in terms of the structure 
functions -Fi and -F2, with introduced when we also consider the Z boson exchange 
contribution: 

[xy^Fi + (1 - y)F2 + y(l - -y)xF3], (2.4) 



with the sum over parton distribution functions now within the structure functions and 
Fi, F2 and F3 receiving contributions from pure 7 and Z exchange and their interference. 

We will now proceed to find the equivalent expression for graviton mediated DIS. The 
advantage of this reformulation of the cross section will be exploited in the next section 
as it facilitates the analytic calculation of the NLO cross section. As with photon DIS 
the starting point is the hadronic tensor. This now has four indices and contains all the 
information about the interaction of the gravitational current with the proton target P 
with spin S: 

W^,^p = 2t? I d^ze'^i^P, S\T^,{z)T^(,mP, S), (2.5) 

with T^i, the energy momentum tensor which describes the coupling of graviton to other 
SM particles. In analogy with the photon case, this hadronic tensor is conserved and can 
be decomposed in term of structure functions. The idea was first introduced in [^], where 
the authors consider a DIS thought-experiment in which the protons are bombarded with 
gravitons instead of leptons. For unpolarised scattering the hadronic tensor is written as: 

3 

i=l 

The tensors A^^^^^ are constructed using the following quantities which are orthogonal 
to the momentum of the graviton, q, so that the conservation condition is automatically 
satisfied: 

P^ = Pf,-{P-q)qf,/q\ (2.7) 

9f,u = 9i^u - q/iQu/q^, (2.8) 
I -2 

id -I) 



P^^P- - 7:r-^9^uP , (2.9) 



with P^j the proton momentum. The relevant tensors are 

= PfiPa9ul3 + PyPa9^l3 + PfiPl39ua + PyPp9^ia 

4 4 — 2_ 

-J—^(.PliPu9af} + PaPpg^^u) + JJ—^P 9f,u9af3, (2-11) 

(3) 2 
^^iual3 = 9^j.a9ul3 +9^il39ua ~ ^ _ ^9^iu9al3- (2-12) 

Each of the tensors above is traceless and orthogonal to the momentum transfer q, two 
conditions that need to be satisfied by the gravitational current. Linear combinations of 
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the structure functions can be obtained by acting on the hadronic tensor with appropriate 
projections. We use as a basis for our projections the following set of tensors: 

Pa = 9ua9ul3 + gpLlSQua, (2.13) 
Pb = PaP^PuPf., (2.14) 
Pc = g,.aP^P^ + guaPf^P^ + g^,pPuPa + gupPf^Pa- (2.15) 

To isolate -Fi, i*2 and -F3 we use the appropriate linear combinations of Pa, Pb and Pc which 
we write down here for simplicity in 4 dimensions, keeping in mind that d-dimensional 
expressions are needed when the NLO corrections are considered within dimensional regu- 
lar isation: 



2x^ 



Pi = -^{Q'Pc - 20x' Pb + 560x^g^PA), (2.16) 
16x2 



P2 = -^{Q'Pc - 2x^Pb + 20x^Q^Pa), (2.17) 
P3 = ^(Q^Pc - x^Pb + 2x*Q*Pa). (2.18) 



The expression for the double differential cross section in Eq. 2.2 is then written in terms 
of the structure functions defined above as 



(fa ^y'^Q^ 



^ (4 - 8y + 52/2 _ y3)Fi - ^(2 - 2y + y^)F2 - 8F3 



dQHx^ 8MI -^-^'^ ' ^ -'-2v- ^ , (2-19) 

where the expressions accompanying the structure functions are obtained by contracting 
the tensors A^^^^p '^ith the leptonic tensor, which only depends on the lepton momentum 
and the momentum transfer. Note that we have rescaled the projection factors so that the 
structure functions are dimensionless. At leading order QCD, taking into account only the 
quark contributions, we have 

= 2x%{x)+q{x)l (2.20) 

= ^[q{x) + q{x)l (2.21) 

while F^ has no quark contribution at LO. 

The corresponding expression for gluon scattering is given by 

^\2-Ay + 2,y^ -y^)g{x), (2.22) 



dQ'^dx 2Ml y2 

which rewritten as above in terms of the structure functions, leads to the following gluon 
contributions to the structure functions: 

Ff = 2x^g{x), (2.23) 
= =^9{x), (2.24) 

F3 = Igix)- (2.25) 

Adding all terms together, we note that at LO Fi and F2 receive contributions from both 
quark and gluon scattering, while F3 only involves a gluonic contribution. 
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2.2 Sum rules 

In this subsection we comment on the existence of sum rules similar to those encountered for 
photon DIS. For photon DIS the Callan-Gross relation F2 = 2xFi is satisfied at LO. This 
sum rule is explained by considering the longitudinal structure function Fl = F2 — 2xFi, 
which is zero at LO as a spin-1/2 quark cannot absorb a longitudinally polarised vector 
boson due to angular momentum conservation. It can be derived by using the polarisation 
vectors of the virtual boson to project out from the hadronic tensor the linear combination 
of Fi and F2 corresponding to the scattering of a longitudinal and transverse photon. At 
NLO the Callan-Gross relation is violated and Fl = 0{as). 

For the graviton, the derivation of a such sum rule at LO is not so straightforward, 
as now both spin-1/2 quarks and spin-1 gluons couple to the graviton. If we now use the 
polarisation tensors of the graviton, formed using polarisation vectors of massive gauge 
bosons we find that the corresponding longitudinal structure function is a linear com- 
bination of Fl', F2 and F^ for which the quark contribution exactly vanishes, that is 
Fj^ oc F^ + IQF2 expected from the same arguments that require F'^ to be zero. Similarly 
by studying angular momentum conservation for the process in the Breit frame we see that 
a quark cannot absorb a graviton with spin component ±2. Therefore the structure func- 
tion corresponding to scattering of a helicity ±2 graviton involves the same combination 
of F^ and F2 ■, so that it receives no quark contribution and also a term proportional to 
F^ . Gluons can absorb any of the five possible polarisations of the graviton so there is no 
sum rule constraining the structure functions based on their gluonic parts. 



3. Graviton exchange NLO results 
3.1 Quark constituent contribution 

Reformulating the cross section in terms of structure functions allows us to extract analytic 
NLO results by simply projecting the structure functions using the NLO hadronic tensor. 
This involves calculating the real and virtual correction diagrams. We start by considering 
initial state quark interactions. The diagrams used to extract the real corrections are 
shown in Fig. |2| and the virtual loop diagrams in Fig. |3|. The calculation of the real 
emission diagrams is first performed at parton level, parameterising the momenta in terms 
of Sudakov variables. For the 0{as) real corrections to the partonic structure functions we 
obtain (with d = A — 2e) 



F\ 



2\e 



Ir 



2-K r(l-e) 



-2x^ 



X 



+ ( 



7 



vr 



)5{l-x) + 



2 - 3x + 



1 / I + 

■ ■ e\ \ — X ^ 
2 - 3x + 3x2 



^5(1 
ln(l — x) 



(1-x). 



1 



Inx 



X 



11 + 0-^^3 

ex e 



1 



1 



X 



2x (1 



-(3 -20x + 20x2) 
(3.1) 
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a) 





Figure 2: Diagrams used to extract tire real di- 
vergences. 





d) 



Figure 3: Loop correction diagrams. 
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+ 



vr 
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74x2 + 20x3) 



^3^ 



2x 

Cpas (4V/Q')^ 



1 



27r r(l - e) 16 



(1 

21 + (1- 



+ 



+ 15x - 12x2 + 3x3 



xj+ 

\2 



1 



Inx 



,(3.2) 



xy 7 - lOx + 3x' 



+41n 



(1 



(2 - 2x + X 



2^^ 



(3.3) 

with /i the renormalisation scale. We note that here we write down only the order Og 
piece of each structure function which has to be added to the LO results presented in 
the previous section. The results above show the expected splitting functions Pq^q and 



q^9 



1 28 1 multiplying the collinear singularities as either of the quark or gluon decay 



products of the initial-state quark splitting can participate in the hard process. We also 
note that the relative prefactor of the splitting functions Pq^q and Pq^g differs between 
Fi and F2. This is directly related to the relative factor of quark and gluon contributions 
to Fi and F2 at LO as shown in Eqs. 2.20| , 2.21, 2.23, 2.24. The absence of the splitting 
function Pq^q from F3 is explained again by the observation that at LO there is no quark 
contribution to F3. 

From the virtual corrections we obtain 



7G 

Iv 



Cpas (4V/Q')'=_2^4 



2 3 
10 



5{1 



(3.4) 



27r r(l - e) 

with diagrams a) and d) of Fig. |3| not contributing within dimensional regularisation. For 
F2 the relative virtual corrections are the same as for Fi, while F3 receives no contribution 
from the loop diagrams. Adding the virtual and real contributions demonstrates the ex- 
pected cancellation of the IR divergences. The divergent terms proportional to the splitting 
functions which correspond to collinear divergences are absorbed into the PDFs. Within 
the MS [29 1 scheme the divergent term, 

1 a, (47r/i2//i2)^ 1 as 

-Pq^qPl = (-7 + 7 - log(47r)) — . 



e27r r(l-e) 



1 

- + 7 
e 



log(47r))^P,^,FfO 



(3.5) 
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where F^^ is the LO partonic structure function at 0{e), is absorbed into the PDF re- 
definition, with the normahsation and factorisation scale here set equal to Q. The form 
of this MS counter term is universal for all processes considered in the remainder of the 
calculation with the splitting function and a possible factor of 2 depending on the process. 

To obtain the hadronic structure functions from the partonic results, we take the 
coefficient functions C(x), defined for example for Fi: 

A = x^gcr(x), (3.6) 

change the argument to z, multiply by q(^^^ / z and integrate over z from x to 1. For 
example the result we obtain for the as part of the NLO hadronic structure function Fi 
involving the quark parton distribution functions is 



vr 



,13 TT . , f dz ,x ^n, |-/ln(l — 2;)\ 2 - 3z + 3z 

(3 - 20z + 20z^) 1 logz 2-3z + 3z 



1 



2z (1-^)+ z 



(3.7) 



with q{x,Q'^) the appropriate NLO parton distribution functions. Here there are no addi- 
tional finite terms originating from considering the LO structure functions at 0(e), as with 
our choice of projections in d-dimensions only and F^ acquire an additional factor of 
(1 -|- e), e.g. FI = — x'^g{x). This has to be consistently taken into account in all the 
contributions to the structure functions that follow in order to extract the correct counter 
term. 



3.2 Gluon constituent contribution 

The same procedure is followed to calculate the gluon contributions to the structure func- 
tions at NLO. The relevant Feynman diagrams for the real corrections are shown in Fig. ^ 
and ^, while the virtual diagrams are shown in Fig. ^, with only the diagrams non- vanishing 
within dimensional regularisation shown. For the real corrections, there exist two sets of 
diagrams with different final states. We note that for the diagrams in Fig. Q we have two 
identical particles in the final states, which implies that we need to divide by two to avoid 
double counting. For all matrix element calculations, when summing over the polarisation 
of external gluons we use 

5;4*(fc)4(^) = -r + - n^^, (3.8) 

in order to include only the physical gluon polarisations. In the results that follow we 
take into account the additional factor of coming from averaging over the initial gluon 
polarisations in d-dimensions. This factor will result in a finite term proportional to the 
splitting function remaining in the coefficient function. 
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Figure 4: Gluon constituent diagrams gG* 
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Figure 5: Gluon constituent diagrams gG* 
qq. 



The results for the partonic structure functions from the diagrams of Fig. I are 
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(3.9) 



2 m2\e 
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25(1 - x) P 
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1 



3e 18e 



(-30 + 2n/)5(l 



^ (l-6)r(l-6) 

1 (38x4 _ 7g^3 _^ 90^2 _ _ 1) _^ 2^ ^^^^ ~ 
(1 — x)+ 6x \ 1 — X y 

^ ^xn ^ ^ (x^-x + l)^ 

+ -5 1 - X - —5 1 - X - 2- ^ logx 

9 3 (1 — x)+ X 



(x^ - X + 1)^ 



(3.10) 



In the results for the real corrections we can identify the Altarelli-Parisi ||2^ splitting 
function: 



l-z 



{i-zy 



+ z{l-z) 



+ i(33 - 2n/)5(l - z) (3.11) 



by rearranging the terms, as 



l-z 



{l-z + z^) 



2\2 



(l-z). 



+ + z{l-z) 



(3.12) 
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Similarly for F2 and F3, 



2 /rt2\e 
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(3.13) 



3r 



_ 3a, (47r/iVQ^)- 1 
27r (1 - e)r(l - e) 8 
1 



2(^(1 - x) 



; 1 + 11.(1 - X) 
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+ ^5(l-x) 



6x 
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(x^ — X + 1) 



-logx 



The corresponding gluon and ghost loop diagram results for Fig. are 



(3.14) 




Figure 6: Initial state gluon virtual correction diagrams. 
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(3.15) 



^2v 



3a, (4^/xVQ2)e 
' 27r (1 - e)r(l - e)y 
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(3.16) 



_ 3a, (4V/Q^)- 1 
27r (1 - e)r(l - e) 8 



3i; 



11 

3^ 



221 

Is 



5(1 -x). 



(3.17) 



Rewriting all the results for the real corrections above, using the expression for the 
splitting function and adding the loop results, we see that we are only left with the collinear 
divergences and the finite terms. Calculating the LO structure function at 0{e) to find 
the appropriate counterterm for the collinear divergence, we obtain the finite coefficient 
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function. Similarly from the diagrams of Fig. I we obtain 



47r (l-e)r(l-e) 
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(3.19) 



(3.20) 



where nf is the number of active flavours. In these expressions, we identify as expected 
the collinear divergence term proportional to the Altarelli-Parisi splitting function Pg^q- 
We note that receives no collinear divergence contribution from the set of diagrams in 
Fig. 1^. This is expected, as at LO there is no quark contribution to i*3. A similar result is 
obtained for the gluonic contribution to for the electromagnetic current. Similarly for 
the fermion loop diagrams of Fig. |6| we have 
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(3.21) 
(3.22) 
(3.23) 



Combining the results proportional to the number of flavours from the loop diagrams 
and the real emission diagrams we verify that the single pole divergences vanish. The 
results for the real and virtual correction diagrams have been obtained using FORM [30] 
and Mathematica, manually performing the Passarino-Veltman [^] reduction for the loop 
diagrams. As a check, a subset of the results were tested using the Mathematica package 



FeynCalc |32| which automatically performs the Passarino-Veltman reduction. 



4. Interference with SM gauge bosons 

As the graviton is a colour singlet, interference effects with other colour singlet gauge 
bosons come into play when we consider the DIS cross section. For the initial state quark 
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interactions we need to consider the interference with the photon and Z boson. At LO the 
G7 interference for positron DIS gives the following contribution to the cross section: 



(4.1) 



(4.2) 



daf^ _ -Avraeg 8 - 12y + 67/^ - 
dQ^ ~ 2M| y ' 

while the GZ interference is 

da^^ ^ X-Ka [dcH^ - I2y + Qy"^ - y"^) - clcljQy - Qy"^ + y^)] 

dQ^ ~ 2Mlsm^29„ (Q2 + m|) y 

where A = ±1 and the sign determines if the interference is constructive or destructive. By 
analogy with the pure photon and pure graviton contributions, the hadronic tensor will be 
of the form 

W^p^ oc J d*z e^'^^P, S\Tp„iz)j40)\P, S), (4.3) 



where is the electromagnetic current. Using the quantities defined in Eq. |2.9| that are 
orthogonal to q, we form a set of three tensors which satisfy the conditions of trace and 
orthogonality, 

y^ppa ^ -papP ^ -ppp^ _ IgP'^p''^ (4.5) 

3 

^^PP- ^ eP^^'I^P'^Po^qp + e^^^^P^Pa^qp. (4.6) 

For the graviton-photon interference only the first two are relevant. These are slightly 
modified for d 7^ 4 to ensure the trace condition is still satisfied, with special care needed 
for the treatment of the 4-dimensional e tensor. When contracted with the leptonic tensor, 
taking into account only the Dirac matrices structure, we obtain the following functions of 
X, y and Q^: 



^(8-12y + 4y2), (4.7) 



Q,n4 

.2 I „,3 



%i-2y' + y'), (4. 



xy^ 

I^(6y-6y2 + y3). (4.9) 

To extract the structure functions we use linear combinations of the following projections: 

PA = PaPpPpand Pb = gppPa + gpaPp, (4.10) 
with Fi and F2 obtained from 

Pi = —^{Q^Pb-20x^Pa) and P2 = -^{Q^Pb - 8x^ Pa). (4.11) 

We can therefore rewrite the differential cross section in the form 

do-T"^ -A7rae„y2 (8 - 12y + 4?/2) (_2y2 + y3) 

dxdQ'^ 2M| y^x^ 
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which at LO gives 



F^^ = x^[q{x) — q{x)\ and F^^ 



(4.13) 



where once more the structure functions are defined to be dimensionless. The quark and 
antiquark contributions to both F^'^ and F2'^ have opposite signs due to the charge factor 
present in the expression. To make the calculation slightly clearer, we can redefine the 
tensor Wi and therefore F^'^ by moving the x^ factor of F^'~' to the tensor Wi which 
leaves the LO partonic F^^ with no explicit factor of x. At the partonic level and again 
using Sudakov parameterisation for the momenta we calculate the diagrams of Fig. 0, where 
the dotted lines denote the final state cut and find 
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Figure 7: Graviton-gauge boson interference diagrams. 
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Similarly for F2 we obtain 
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We note that apart from an overall factor, the two expressions differ only in the finite 
terms. Calculating the interference terms between a graviton exchange diagram at 0{as) 
and a LO order photon exchange diagram and vice versa we obtain the following expression 
for the virtual corrections: 



^7G _ Cpas (47r/LiVQ 
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(4.16) 



The constant term originates from combining the constant term contribution from the 
virtual loop in photon exchange (—8) and graviton exchange (—10). The result for ^2'^ is 
identical. Using 



(l + x^ 



1 



l + 3 
737 -2'^^'-")' 



(4.17) 



(l-x)_ 

we notice that as expected the double and single pole divergences exactly cancel when 
we add the virtual and real contributions, and that we obtain the expected form of the 
collinear divergence: 



27r r(l-e) e^'^'- 



The relevant coefficient function is 
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(4.19) 
(4.20) 
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We note that there is a second set of interference diagrams with initial state gluons shown 
in Fig. ^. We calculate the matrix element squared for this contribution; however, the 
result is odd in cosO* where 9* is the angle between the outgoing quark and the initial state 
gluon in the gluon-graviton( /photon) centre-of-mass frame. Therefore when we perform 
the two-body phase space integration we get no contribution to the structure functions F^^ 



and -^2^*^. A similar effect arises in [21| where the authors note that there is no interference 



between the graviton and SM contribution to Drell-Yan. This does not occur at the matrix 
element squared level but only when the phase space integration is performed. 
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Figure 8: Interference diagrams with initial state gluon. 

For the Z-graviton interference the odd parity structure function comes into play. The 
projection used to extract it is 



Ax 



(4.23) 



with the other two projections remaining the same. The differential cross section is now 
written as 



da 



ZG 



where at LO 



TraA 



2, ,2 



dxdQ^ 2M|sin22e„ {Q^ + M|) 
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(4.24) 



pGZ 
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x^cl[q{x) - q{x)], 



[q{x) -q{x)], 



X 

3 - -^[Qix) + Qi^)]- 



GZ 



(4.25) 
(4.26) 

(4.27) 



As the first two structure functions are independent of 7^ and given the universality of the 
0{as) QCD corrections for spin-1 gauge bosons, we expect the NLO structure functions 
^ and ^ to receive the same relative corrections as and F2^ respectively. For 
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the NLO contribution to F^^ we need to consistently treat 7^ in d dimensions. Following 
the prescription in [33|, we replace the axial current by 



(4.28) 



and then perform the trace with the gamma matrices in d-dimensions and keep the product 
of the two e factors (one from the projection and one from the matrix element) outside 
the trace. The product of epsilon tensors gives a determinant of Kronecker deltas which 
are then treated as d-dimensional objects. Following consistently this prescription for the 
virtual and real corrections, we obtain for the real corrections to the partonic structure 
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(4.29) 



while from combining the loop corrections to Z and graviton vertices we get for the virtual 
corrections 



Fi 



8tt 



r(i 



l-i-4 

e2 e 



5{l-x), 



(4.30) 



and therefore the divergences cancel as expected. 



5. Numerical results 



The analytic results for the cross sections are implemented into a program to numerically 
integrate over the momentum fraction and the DIS variable x. The integration limits for 
X are normally to 1 with the lower limit modified to Q^/O.Os because experimentally 



at HERA the maximum value for y is 0.9 |34|. The proton energy is 920 GeV and the 
positron energy is 26.7 GeV and we calculate the differential cross section in the range of 
experimentally measured values of Q^. 

We first show the effect of the NLO corrections to pure graviton exchange, splitting 
the cross section into contributions from intitial state quarks and gluons. The calculation 



is performed using the NLO PDF set MSTW2008NLO [|5[, with the factorisation and 
renormalisation scales set to Q. In Fig. ^ we show the LO and NLO results for the differ- 
ential cross section da/dQ"^ for Ms = 1 TeV and the HERA positron and proton energies. 
This result is proportional to AfJT^, so one can appropriately rescale for other values of 
Ms- In Fig. 10 we decompose the results to gluon and quark initial state contributions 



while in Fig. |ll| we show the relevant /c— factors where k is defined as the ratio of the 
NLO cross section calculated with NLO PDFs to the LO cross section calculated with LO 
PDFs. We notice that the NLO corrections reduce the total cross section by up to 10%, 
with the reduction being more significant for the gluon processes. This is due to the fact 
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Figure 9: Differential cross section for graviton exchange for Ms — 1 TcV. 



that the 0{as) corrections are dominated in both cases by large and negative contributions 
proportional to 5{\ — x). For both LO and NLO cross sections, the dominant contribution 
comes from quark scattering, as we are considering the cross section at high values of the 
momentum fraction x where the gluon PDF suppresses the gluon contribution. 
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Figure 10: Relative contributions to the cross 
section from initial state quarks and gluons. 
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Figure 11: Relevant fc— factors for the contribu- 
tions from initial state quark and gluons. 



The advantage of the NLO calculation is that the result for the cross section is less 
dependent on the factorisation and renormalisation scales, and therefore the theoretical 
uncertainty is reduced. The dependence of the results on the scale choice is calculated 
by reinstating the appropriate \og{Q'^ / ^'^) factors in the coefficient functions. This can be 



seen in Figs. |12| and |13| where we show the dependence of the LO and NLO cross section on 
the choice of factorisation and renormalisation scale which here are chosen to be equal and 
varied as a multiple of the momentum transfer Q. We consider two values of to show 
the differential cross section for M5 = 1 TeV, = 250 GeV^ and 12000 GeV^. Similar 
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results are obtained for all values of Q^, and therefore the total cross section follows the 
same behaviour in terms of the scale dependence. 
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Figure 12: Scale dependence of the differential 
cross section da/dQ^ at = 250 GeV^. 



Figure 13: Scale dependence of the differential 
cross section dcr/dQ'^ at = 1200O GeV^ 



The effect of the NLO QCD corrections is also considered for the interference terms. 



The /c— factor in this case is found to be very close to one as shown in Fig. 14 for 
Ms = 1 TeV. This is due to the fact that the pure signal contribution receives a neg- 
ative contribution from NLO corrections while the SM prediction is increased by a similar 
relative amount when we consider order a, corrections. 
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Figure 14: Differential cross section from interference terms for Ms — 1 TeV, and the correspond- 
ing fc— factor. 



In Fig. 15 we show the ratio of the total differential cross section da/dQ^ for the NLO 
new physics prediction to the NLO SM prediction. The data shown here are taken from 
1 34]. The experimental errors on the data points are the total errors given in column 7 of 
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Table 6 in [34|. We note that for Ms=l TeV we would only observe the interference effects 
which have opposite sign for A = ±1. At sufficiently high the pure graviton contribution 
takes over, giving a positive deviation from the standard model prediction. However the 
cross section falls rapidly with and the values fall beyond measurable event rates. 
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Figure 15: Ratio of differential cross sections for Ms = 1 TeV. 



Given the recent LHC results from dijet measurements [0, 12] and the subsequent 
analysis of data which excludes fundamental scales up to ~ 3 TeV |13|, the graviton 
contribution to DIS at HERA energies could not be detected. While the reach of HERA 
experiments is limited by the low centre-of-mass energy, the situation is more promising for 



the proposed LHeC. In Fig. 16 we show the cross section for collisions of 7 TeV protons and 
140 GeV electrons, keeping the same maximum y value as at HERA for comparison with 
Fig. ^ and the relevant fc— factor. We note that the A:— factor does not vary significantly 



between the two different centre-of-mass energies. In Fig. 17, we show the deviation from 
the SM prediction for Ms = 1 TeV for comparison with the prediction in Fig. 15 for HERA 
and to show the extended reach of the experiment also for Ms = 3 and 4 TeV. We use 
different colours to plot the results for different scales and solid and dotted lines to show the 
effect of changing the sign of the interference terms. We notice that even for a fundamental 
scale as large as Ms = 4 TeV, we can get a measurable deviation from the SM prediction, 
with the cross-section values now much higher leading to measurable event rates. Note 
that as the sign of the 7G interference term is different the picture is slightly different 



from Fig. 15, as the first deviation (lowest Q^) from the SM prediction comes from the 
interference term. 



6. Conclusions 

In this paper we have studied the graviton contribution to DIS, defining appropriate struc- 
ture functions and obtaining analytic results for the order differential cross section. We 
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Figure 16: Differential cross section for tlie 
LHeC at LO and NLO for Mg = 1 TeV and the 
relevant fc— factor. 



Figure 17: Ratio of the NLO new physics and 
SM differential cross sections at LHeC energies for 
Ms = 1,3,4 TeV. 



calculate the relevant /c— factors and find that NLO QCD corrections decrease the cross 
section by up to 10% with quark scattering contributing most of the cross section in the 
relevant region of high Q^. We also show that, as expected, the NLO results are less de- 
pendent on the choice of factorisation and renormalisation scales and therefore decrease 
the theoretical uncertainty. We conclude that NLO corrections would only slightly mod- 
ify the limits set on the fundamental scale, compared with LO calculations, as the limits 
depend mostly on the interference of the SM and graviton contributions. Given currents 
constraints set recently by LHC results, the graviton contribution to DIS is not within the 
HERA reach. However, in the context of the proposed LHeC experiment, given the signifi- 
cantly higher center-of-mass energy the graviton contribution to DIS could be observed for 
scales as high as 4 TeV. We have calculated the cross section and the expected deviation 
from the SM prediction at NLO with the calculated fc— factor showing a similar behaviour 
as for the lower HERA energies. 
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